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Abstract With the aim of understanding the role of massive outflows in high-mass star formation, we mapped 
in the ^^CO J = 2 — 1 transition 26 high-mass star-forming regions at very early stages of their evolution. At 
a spatial resolution of 11” bipolar molecular outflows were found in 21 of them. The other five sources show 
confusing morphology but strong line wings. This high detection rate of bipolar structure proves that outflows 
common in low-mass sources are also ubiquitous phenomena in the formation process of massive stars. The flows 
are large, very massive and energetic, and the data indicate stronger collimation than previously thought. The 
dynamical timescales of the flows correspond well to the free-fall timescales of the associated cores. Comparing 
with correlations known for low-mass flows, we And continuity up to the high-mass regime suggesting similar 
flow-formation scenarios for all masses and luminosities. Accretion rate estimates in the 10^ Lq range are around 
10 “^ M 0 yr~^, higher than required for low-mass star formation, but consistent with high-mass star formation 
scenarios. Additionally, we And a tight correlation between the outflow mass and the core mass over many orders 
of magnitude. The strong correlation between those two quantities implies that the product of the accretion effi¬ 
ciency /acc = Afacc/(ALcore/tff) and fr (the ratio between jet mass loss rate and accretion rate), which equals the 
ratio between jet and core mass (/acc/r = Mjet/Mcore), is roughly constant for all core masses. This again indicates 
that the flow-formation processes are similar over a large range of masses. Additionally, we estimate median fr and 
/acc values of approximately 0.2 and 0.01, respectively, which is consistent with current jet-entrainment models. 
To summarize, the analysis of the bipolar outflow data strongly supports theories which explain massive star for¬ 
mation by scaled up, but otherwise similar physical processes - mainly accretion - to their low-mass counterparts. 

Key words. Molecular data - Turbulence - Stars: early type - Stars: formation - Stars: mass-loss - ISM: jets and 
outflows 


1. Introduction 


Molecular outflows are a well known phenomenon in sites 
of low-mass star formation, and they play an important 
role in transporting angular momentum away from the 
forming star. The observational database on low-mass out- 
flows has increased tremendously over the last decade, giv- 
ing rise to different formation scenarios. While it is still 
unclear how the outflow is accelerated near the proto¬ 
star and/or disk, it is now widely believed that low- 
mass molecular flows are momentum driven by highly 
collimated jets, which entrain the surrounding molecular 
gas. For recent reviews on this topic see, e.g., [Richer et 


al. (2000 ); Bachiller &: Tafalla (2000| ); $hu et al. (200Cj) ; 
Konigl fc Pudritz (200C ). 

The situation in the high-mass star formation regime 
is less clear, because, due to the rarity of these objects and 
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the typically larger distances (a few kpc) of sites of high- 
mass star formation, the spatial resolution has been lack¬ 
ing to resolve the outflows and their driving sources prop¬ 
erly. Recent systematic molecular line studies of high-mass 
outflow s h ave been carried out by Shepherd fc Churchwelj 
(1996a , ^); Henning et al. (2000 ); Zhang et al. (2001 ); 
Ridge fc Moore (2001 ). All of them And that ~ 90% of 
the observed massive star-forming regions are associated 
with high v elocity gas. Mapping a subsample of 10 sources 
(out of 94) Shepherd &: Churchwell (1996b) And bipolar 
morphology in 5 of them and Zhang et al. (2001 ) found 
spatial outflow structures in 39 of 69 sources. Thus, both 
studies indicate bipolar structures in at least 50% of the 
observed sources. 


Including more massive outflow sources from the lit¬ 
erature, the overall picture has emerged, that outflows 
are ubiquitous phenomena in massive star formation, that 
they are very massive (up to hundreds of solar masses in 
the flows), very energetic (~ 10'^® erg) but seemingly not 
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very collimated (i.e. collimation factors - the length of the 
flow divided by its width - between 1 and 10 for low-mass 
sources versus collimation factors between 1 and 1.8 for 


high-mass sources, [Richer et al. 200(]| ; |Churchwell 2000b| ; 
Ridge fc Moore 2001 ). The high masses in the flows as 
well as the low collimation factors, which are difficult to 
explain in current jet-entrainment scenarios, have given 
rise to new ideas of outflow formation including the pos¬ 
sibility that massive outflows consist of accelerated gas 
which has been deflected by the young accreting proto¬ 


star, rather than swept-up ambient material (Churchwell, 
2000b| ^ evine et al. (199S| ) propose a slightly different 
scenario in which massive stars produce collimated jets 
only in their earliest phases, and with increasing lumi¬ 
nosity of the central objects the jets may be replaced by 
wide-angle winds. 

To test the proposed scenarios it is necessary to in¬ 
crease the number of outflow sources using a homogenous 
sample with adequate angular resolution. Therefore, we 
mapped 26 promising candidates from a larger sample of 
massive star formation sites in the J = 2 ^ 1 line of CO 
with the 30 m telescope at a spatial resolution of 11". This 
work is part of a long term project to search and investi¬ 
gate 69 massive and rather isolated star formation regions 
in an evolutionary stage prior to building significant ul¬ 
tracompact Hii regions. An introduction to the sa mple 
and first results are g iven in ^ridharan et al. (2002 ) and 
Beuther et al. (2002a ). 

In the following, we present the observed dataset and 
give a detailed analysis of the outflow characteristics. The 
derived parameters are compared with those from low- 
mass as well as other high-mass outflows, and we discuss 
the implications for star formation in such regions. 


2.2. IRAM 30 m observations 

The IRAM 30 m telescope on Pico Veleta near Granada 
(Spain) was used to map 26 sources in the J = 2 — 1 
transition of in April and November 1999 and in 

April 2000. Maps with an average extent of 100" x 100" 
were observed in the on-the-fly mode, where the telescope 
moves continously across the source and dumps the data 
in a well chosen grid. At 230.5 GHz the resolution of the 
30 m telescope is approximately 11". Beuther et al. (2000 ) 
have shown that Nyquist-sampling is sufficient in the on- 
the-fly mode as well, because the beam is smeared out by 
only about 4%. Therefore we chose as a sampling interval 
4" (Nyquist interval ^ ^ 4.5"). The dump time 

was 2 sec per position and we scanned all maps twice, 
in most cases in perpendicular directions to reduce scan¬ 
ning effects. Experience showed that the whole system 
(weather+technical system, Tsys ~ 250 K) is stable for 
around 5 minutes. Therefore our ON-OFF cycles never 
exceeded that time. OFF positions were chosen based on 
the Stony Brook and the CfA-survey of the galactic plane 
( ganders et ah, 1986 ; Dame et ah, 2001) and checked to 
be emission free. The frequency resolution was 0.1 km/s 
and the beam efficiency 0.41. 

The data were reduced with CLASS and GRAPHIC 
of the GILDAS software package by IRAM and the 
Observatoire de Grenoble. To improve the signal to noise 
ratio, the ^^CO data are smoothed to a velocity resolution 
of 1 km/s, sufficient to sample the broad CO lines. 


2. Observations 

2.1. Source Selection 

The 26 sources observed in this study were chosen from a 
sample of 69 high-mass proto-stellar candidates based on 
the CO 2-1 li ne wings seen by thes e authors in pointed 
observations ( Sridharan et ah, 2002| ). We do not believe 
that there was a particular bias in this selection, ^ridharan 


et al. (20021 ) found wings in 85 percent of their sample 
and argue on the basis of the statistical distribution of 
inclination angles that essentially all of these sources are 
associated with outflows. As also argued by ^ridharan et 


al. (2002|), the majority of this sample is likely to be young 


and 50 percent have no associated 9 GHz radio continuum 
flux down to a limit of 1 mJy. This suggests they have not 
had time to form a substantial Hii region. While the low 
radio continuum flux may have a variety of causes (e.g., 
high 9 GHz optical depth or dust mixed with ionized gas), 
we consider the sample of sources observed by us to be an 
arbitrarily chosen set of very young regions of high-mass 
star formation. 
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Figure 1. Massive molecular outflows: the solid lines show the blue wings and the dashed lines the red wings of CO 
2-1 emission. The contours are chosen to highlight the most prominent features in each source, usually between 30% 
and 90% (steps of 10%) of the peak integrated wing intensity. The dotted lines in 18566+0408 represent an additional 
feature at “blue” velocities [58,60] km/s. The grey scale p resents the correspond ing 1.2 mm continuum maps from 5% 
and 95% by steps of 10% of the pea k flux as outlined in [Beuther et al. (2002a ). The axis show offsets in arcsec from 
the absolute IRT5-positions given in Sridharan et al. (2002|). Detailed morphological descriptions are given in Sect. 3.1 
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Figure 1. continued. The red wing for the southern flow in 20293+3952 is missing due to instrumental problems in 
that region, but simultaneously observed SiO data (not presented here) confirm this southern flow. 
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Figure 2. For sources without clear bipolar outflow structure we show the summed 2-1 spectra (Tmb [K] versus 

velocity [km/s]) corresponding to the regions outlined by the 1.2 mm continuum maps in grey scale. The continuum 
images are presented as in Fig. . 
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3. Observational results 

3.1. Morphologies 

Fig. H presents the outflow sources [CO 2-1 blue (solid 
lines) and red wing emission (dashed lines) overlayed on 
the grey scale 1.2 mm continuum maps (Beuther et al 


2002a)]. From the 26 sources observed, 21 show either 


bipolar morphology or wide non-Gaussian line wings such 
that we suspect the outflow is along the line of sight 
(23139+5939) and not spatially resolvable with the 11" 
resolution of our observations. The 5 remaining sources 
(18089-1732,18308-0841,19220+1432,19413+2332 and 
20205 + 3948) also show line wings, but the maps are very 
chaotic, so that we could not define any o utflow structure. 
All the sources are in the galactic plane (Sridharan et ah. 


2002 ) and unrelated spiral arm emission is partly confus¬ 
ing the spectra. Fig. || presents for “confused” sources the 
sum of all s pectra in the regions o utlined by the 1.2 mm 
dust maps ( Beuther et ah, 2Q02a ). We cannot tell if the 
latter are a superposition of many flows, or if they are 
more isotropic flows, similar to that observed in Orion 
( Rodriguez-Franco et ah, 199S ; Schultz et ah, 1999 ). This 
can only be investigated with higher angular resolution 
in the future but we note that most of our sa mple is less 
luminous than Orion IRC2 (L ^ 2 x 10® Lq, 


Churchwell 


2000b|) . Within the range of luminosities which we sample 


(roughly 10® ® to 10® Lq), we And that bipolar features 
are omnipresent. We also do not observe a statistical dif¬ 
ference in finding bipolar outflow structures for sources 
with or without detected cm emission suggesting perhaps 
that the arrival of proto-stars on the ZAMS is not critical 
for the outflow. 


the mm peak (18566+0408, 19035+0641, 19217+1651, 
19411+2305, 22570+5912). The reasons for such offsets 
are not clear. For 18566+0408 and 19035+0641, pointing 
problems may be the explanation since simultaneously ob¬ 
served 1-0 maps (not published so far) peak ap¬ 

proximately at the center of the flows and offset relative 
to the mm cores. In the cases of the more chaotic looking 
flows 19411+2305 and 22570+5912, this seems unlikely 
and the structure may be due to several low-mass flows, 
which are not at the center of the mm core. 

19217+1651 on the other hand is not chaotic but a nice 
examples of a bipolar outflow which makes a superposi¬ 
tion of several low-mass flows very unlikely. Additionally, 
the offset between the flow center and the single dish out¬ 
flow map is real because H^®CO+ 1-0 as well as thermal 
CH 3 OH (both observed simultaneously with the outflow) 
peak at the mm core. Using a theoretically derived cor¬ 
relation between the moment um of the outflow an d the 
proto-stellar mass (Eq. 18 in [Tan fc McKee (2002 )), we 
can infer a mass for the proto-star powering the outflow. 
In this way, we conclude that the proto-star has a mass of 
at least 10 M©. This suggests that high-mass proto-stars 
can emerge from the core where they were formed. The 
projected offset between the ffow center and the mm peak 
is of order 0.3 pc which implies a relatively high peculiar 
velocity of the proto-star relative to the core (~ 4 km/s) 
using the ffow timescale given in Table ||. This is obviously 
an object of considerable interest and interferometric ob¬ 
servations with the Plateau de Bure Interferometer are 
currently being conducted. 

The flows are rather large with an average size of 
~ 0.8 pc (see Table ||). The collimation factors derived 
for this sample are higher than previously claimed. Richer 


vious sing le pointing observation s conducted at the au m 
telescope ( Sridharan et ah, 2002| ). Therefore we expected 
outffows in all sources, but the large fraction with bipolar 
structures (~ 80%) is rather surprising. Previous mapping 
studies of massive molecular outffows found bipolar struc¬ 


All the mapped sources showed wing emission in pre et al. (2000) in their review note that so far no massive 


tures only in approximately 50% of the sources ( [Shepherd 
& Churchwell, 1996a; Zhang et ah, 2001), which seemingly 
suggested that bipolarity is less probable in high-mass 
than in low-mass star formation (Churchwell, 2000b). The 
high detection rate of bipolar outffows in our study is most 
likely due to the improved spatial resolution of our data 


(11") (homparecl to previous statistical work on massive 
outffovre (compare radius r m Ta ble [I|to 60 " resolution m 


Shepherd fc Churchwell (1996b ) and 30" in Zhang et al. 


(2OOII )). 


Every outflow is associated with massive mm cores 
(grey-scale in Figure |^). It is not evident that the size 
of the flows depends on the size of the cores: about 
50% of the flows have an extent larger than the mapped 
core sizes, whereas the other half of the flows are 
smaller. In most cases, the ffows are centered on the 
mm peak, which most likely harbors the most massive 
and youngest proto-stars. But there are a few remark¬ 
able exceptions, where the ffows are centered offset from 


ffow with a collimation factor fc larger than 1.8 has been 
observed. The last column of Table presents the col¬ 
limation factors we derive for our sample (ffows with a 
question mark are along the line of sight, which makes 
a determination of fc difficult). These collimation factors 
are lower limits due to angular resolution and projection 
effects, but we nevertheless obtain a mean fc of 2 . 1 , larger 
than most measurements of high-mass ffows to date. The 
largest fc we ffnd in the single dish data is 3, but we know 
from interferometric observations of selected s ources that 
fc can be even higher: fc ^ 10 for 05358+3543 ( [Beuther 
ah, 20Q2b ) and /c ~ 5 for 23033+5951 (Wyrowski et ah, 
in prep.). It is interesting to compare the observed col¬ 
limation factors with theoretical expectations: a ffow at 
5 kpc (the mean distance of this sample) extended 1 pc in 
length and 0.1 pc in width results in an intrinsic fc = 10 . 
Inclination effects to first order just reduce the length of 
the ffow but not the widths; thus the length has to be 
corrected on the average by cos(57°) (mean inclination 
angle). After converting the linear scales into arcseconds 
and convolving with an 11" beam - as for our CO data 
- the length of the ffow is ^ 25" and the width ~ 12", 
resulting in an observable fc = 2 . 1 , which is similar to the 
mean value we derive for our sample. Thus, the observed 
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collimation factor is a strong function of the available an¬ 
gular resolution and high-mass flows may in general be as 
well collimated as flows from low-mass proto-stars (though 


see the discu ssions of Devine et al. (1999 ) and Reipurth fc 
Bally ( 2001 )). Even this new dataset with 11" resolution 
underestimates collimation factors; higher resolution with 
interferometers is needed for that. 

In 4 fields (18151-1208, 19410+2336, 20293+3952, 
20343+4129) we see at least two separated flows. 
Interferometric data for two of the flows (05358+3543 
Beuther et al. 20021; 23033+5951 Wyrowski et ah, in 


nation angles of t he flows ar e also often unknown. [Cabrit 


(rrmb P^C02-l)7 meaub & meaur) were used as inputs 
to the flow calculations. Additional input parameters are 
the radius r of the flow from the projected center, and 
the maximum velocity separation of the red and the blue 


prep.) reveal that those single dish flows split up into at 
least two outflows, but the data also show that most of 
the single dish flux (> 50% for 05358+3543 and > 80% 
for 23033+5951) is produced by one major flow with minor 
contributions from the second flow. There are likely to be 
m ultiple lower-mass flo ws, and recent theoretical studies 
by Tan fc McKee (2002 ) consider massive outflows to be a 
superposition of outflows from all proto-stars of each clus¬ 
ter evolving simultaneously. This scenario results in less 
collimated flows and may explain some of our more chaotic 
sources, but it does not seem to be applicable in general. 
Thus, contributions from lower-mass flows are possible, 
but the high degree of collimation and bipolarity we And 
on average in addition to that found in a small number of 
sources with interferometric observations strongly suggest 
that our data are dominated by the most massive outflow 
of each cluster. 


3.2. Characteristics of the outflows 

The determination of outflow parameters is subject to a 
large number of errors. To first order, it is difficult to sepa¬ 
rate the outflowing gas from the ambient gas, and the incli- 


wings nniaxr & 'I'maxb- All the input parameters are pre¬ 
sented in Table An additional uncertainty is caused by 
distance ambiguity in which cases the outflow parameters 
are calculated for the near and the far distance (Table ||) . 

Opacity corrected H 2 column densities N\, and in 
both outflow lobes can be calculated by assuming a con¬ 
stant ^ ^CQ/^^CQ 2-1 line wing ratio throughout t he out¬ 
flows (Cabrit & Bertout, 1990). Choi et al. (1993| ) found 
an average ^^CO/^^CO 2-1 line wing ratio around 0.1 in 
7 massive star-forming regions, which we adopt for our 
sample as well. Levreault (198^ ) found similar values in 
low-mass outflows. The derived flow parameters are the 
masses Mb, Mr in the blue and red outflow lobes and 
the total mass Mout, the momentum p, the energy E, the 
size, the characteristical time scale t (radius of the flow r 
divided by the flow velocity), the mass entrainment rate 
of the molecular outflow Mout (as opposed to the mass 
loss rate of the jet Mjet), the mechanical force Fm and 
the mechanical luminosity L m for each flow. For details 
see Cabrit fc Bertout (199( ). The used Eqs. are (assum- 


ing H2/^^C0 = 89 X 10"^ ( Cabrit fc Bertout, 1992| ) and 
Tex = 30 K; ruHa is the mass of the H 2 molecule): 


fc Bertout (1986| , |l990| , [1992[) performed a series of studies 
on outflows from low-mass sou rces and the approach w e 
are following here is outlined in Cabrit fc Bertout (1990|) . 
The velocity range due to the outflow is determined in 
two ways. One is to define the flow by the line wings in 
the spectra, the other approach is to map each channel 
and decide from the spatial separation of different chan¬ 
nels, which belong to an outflow and which correspond 
to the ambient core emission. In some sources both crite¬ 
ria were suitable, in other sources only one of them, e.g., 
the 0535 8+3543 flow is know n to be near the plane of 
the sky ( Beuther et ah, 2002b ) causing us to use the spa¬ 
tial separation. A counter example is 23139+5939, which 
shows extremely strong wing emission but nearly no lobe 
separation probably due to being rather along the line of 
sight (Fig. 0 )- Then we mapped the chosen wings, and 
the size of the flow (transformed to linear scales, sizeb 
fc sizer) and the mean value of integrated wing emission 


N 

Mout 

P 

E 

t 

Mout 

Fur 
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= ( 


H 2 


//T„rb('"C02 - 1) 


y^COj V/r„rb(i2C02-l)y 47r3M2^2 


X e 


-16.6/Te, 


^ Tr„b(''C02 - 1) 

(A^b X sizeb + Nr X sizer) mua 

Mb ^ 'l^maxb + Afr X 'Umaxr 


-Mb X 
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+ 2'^'' ^ ^maxr 


(V 


'maxb 
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M, 
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P 

t 

E 
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( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


In the case of 18151 — 1208, we just determined the char¬ 
acteristics for the eastern flow, because the western one 
is confused by ambien t emission. The derived quantities 
are shown in Table ||. Cabrit fc Bertout (199C| ) state the 
mass determinations to be correct within approximately a 
factor 2, while kinematic parameters should be approx¬ 
imately correct within a f actor 10. It is clear however 
from the Choi et al. (1991^) study that the high opacity 


in ^^CO transitions can lead to large errors for individ¬ 
ual objects and it would be useful to obtain data 

for these sources. While several of the derived parameters 
are discussed in detail in we want to stress that all of 
the observed sources show massive and energetic molec¬ 


ular outflows compared to low-mass flows (Richer et al. 


200 c). 


Another important result from our study concerns the 
derived flow timescales t. These are likely to represent a 


lower limit to the true age of the flows (see Parker et al 
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Table 1. Parameters f or outflow calculations : The number of flows per source the infrared luminosity L derived 
from the HIRES-daXa. ( Sridharan et al., 2002 ) (if the kinematic ambiguity is unresolved both values - near and far - 
are given), the velocity range of the wings An, the mean of the integrated wing emission, the maximum projected 
velocity and the radius r of the flow. The last column presents the collimation factor fc 


source 

# 

L 

[log(L 0 )] 

Aub 

[km/s] 

AVt 

[km/s] 

meanb 
[km/s K] 

meaur 
[km/s K] 

Umaxb 

[km/s] 

"^^maxr 

[km/s] 

r 

n 

/c 

05358-73543 

1 

3.8 

(-32,-21) 

(-12,-4) 

28.5 

29.7 

14.4 

13.6 

60 

2.0 

18151-1208 

1 

4.3 

(25,30) 

(38,42) 

20.2 

12.2 

5 

10 

45 

2.1 

18182-1433 

1 

4.3/5.1 

(53,56) 

(62,70) 

26.5 

26.1 

6.1 

12 

25 

1.7 

18264-1152 

1 

4.0/5.1 

(28,37) 

(52,63) 

38.4 

46.6 

16 

18 

20 

? 

18345-0641 

1 

4.6 

(82,90) 

(102,106) 

28.6 

16.8 

14 

10 

40 

1.5 

18470-0044 

1 

4.9 

(86,91) 

(102,106) 

24.3 

9.0 

10.5 

9.5 

38 

2.0 

18566-70408 

1 

4.8 

(68,77) 

(96,102) 

15.2 

9.2 

17 

17 

40 

2.4 

19012-70536 

1 

4.2/4.7 

(55,62) 

(70,75) 

39.6 

21.6 

10 

10 

30 

? 

19035-70641 

1 

3.9 

(21,27) 

(38,43) 

12.0 

4.2 

11 

11 

30 

? 

19217-71651 

1 

4.9 

(-16,-3) 

(9,21) 

31.4 

26.2 

19.5 

17.5 

28 

2.8 

19266-71745 

1 

1.7/4.7 

(-3,0) 

(10,15) 

15.4 

5.4 

8 

10 

35 

? 

19410-72336 

1 

4.0/5.0 

(1,17) 

(26,45) 

35.3 

64.0 

21.4 

22.5 

43 

2.0 

19410-72336 

2 

4.0/5.0 

(1,17) 

(26,45) 

36.0 

52.0 

21.4 

22.5 

45 

2.1 

19411-72306 

1 

3.7/4.3 

(15,25) 

(35,43) 

32.4 

27.0 

14 

14 

35 

2.3 

20216-74107 

1 

3.3 

(-9,-4) 

( 1 , 6 ) 

30.8 

28.8 

6.6 

8.2 

32 

? 

20293-73952 

1 

3.4/3.8 

(-16,4) 

(12,25) 

25.6 

61.6 

26 

30 

32 

2.7 

20293-73952 

2 

3.4/3.8 

(-16,4) 

(12,25) 

18.2 

3.4 

26 

30 

15 

? 

20343-74129 

1 

3.5 

(8,9) 

(13,15) 

24.6 

33.4 

3.5 

4.5 

23 

1.0 

20343-74129 

2 

3.5 

(8,9) 

(13,15) 

9.2 

12.0 

3.5 

4.5 

25 

1.7 

22134-75834 

1 

4.1 

(-30,-23) 

(-14,-4) 

10.6 

34.4 

11.7 

14.3 

37 

? 

22570-75912 

1 

4.7 

(-56,-50) 

(-42, -32) 

20.8 

17.8 

10 

14 

60 

? 

23033-75951 

1 

4.0 

(-67,-60 ) 

(-47,-33) 

21.8 

34.2 

14.0 

19.4 

50 

3.0 

23139-75939 

1 

4.4 

(-59,-47) 

(-39,-31) 

43.0 

23.6 

11.1 

13.2 

29 

? 

23151-75912 

1 

5.0 

(-89, -65) 

(-47, -35) 

53.8 

11.8 

34 

19 

20 

1.2 


1991) and hence also to the time over which the embedded 
proto-stars responsible for the flows have been accreting 
from their surroundings. One can compare these estimates 
for the outflows with the free-fall times (ts = 
derived from the mean gas density p in the cores from 
which presumably the accretion takes place, p is derived 


from t he dust continuum maps as discussed by Beuther 
et al. (2002a ). The free-fall times are an estimate of 
the timescale for dynamical evolution of the core and in 
fact theoretical studies of core evolution suggest that the 
timescale for star formation in a core is of order a few free- 
fall times (e.g.. Tan fc McKee 20^ ). Our comparison of t 
and Iff is shown in Fig. |^. We see that flow timescales and 
free-fall timescales are of the same order of magnitude and 
range for our sample between 2 x lO'^ and 2 x 10® years with 
an average dispersion about a factor 2 (for the sources with 
known distance). The rough equality between these inde¬ 
pendently derived quantities is surprising and it will be 
important to check whether this also holds for other sam¬ 
ples. We tentatively conclude that for the present sample, 
flow ages are good estimates of proto-star lifetimes and 
will assume this in the following discussion. 



Figure 3. Dynamical outflow times t versus free-fall times 
tff. The circles show sources with known distance whereas 
triangles represent sources with unresolved distance am¬ 
biguity (open triangle: near distance, fllled triangle: far 
distance). 


4. Discussion 


4.1. High-mass versus low-mass outflows 

We wish first of all to compare our results to systematic 
studies conducted in the low-mass regime. The analysis 
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Table 2. Outflow results: distances D [kpc], number of outflows column densities N [logio(cm~^)], masses 
Mb (blue), Mr (red) and Mout (Mout = Mb + Mr) [Mq], momentum p [Mq km/s], energy E [lO^^erg], size [pc], 
time t [lO^yr], mass entrainment rate M^out [10~'^M©/yr], mechanical force [lO“^M 0 km/s/yr] and mechanical 
luminsoity Lm [i©] 


source 

D 

# 

Nb 

Nr 

Mb 

Mr 

Mout 

P 

E 

size 

t 

Mout 

Fr. 

Lra 

05358-t3543 

1.8 

1 

21.1 

21.2 

11 

10 

21 

288 

4.0 

0.52 

3.7 

5.6 

7.9 

9.1 

18151-1208 

3.0 

1 

21.0 

20.8 

8 

4 

12 

106 

0.9 

0.65 

7.1 

1.7 

1.5 

1.1 

18182-1433 

11.7 

1 

21.1 

21.1 

130 

73 

203 

1665 

15.0 

1.42 

15.3 

13.0 

11.0 

8.2 

18182-1433 

4.5 

1 

21.1 

21.1 

19 

11 

30 

246 

2.3 

0.55 

8.9 

3.4 

2.8 

2.1 

18264-1152 

12.5 

1 

21.3 

21.4 

234 

171 

405 

6818 

110 

1.21 

7.0 

58.0 

98.0 

135.0 

18264-1152 

3.5 

1 

21.3 

21.4 

18 

13 

31 

534 

9.0 

0.34 

3.7 

8.6 

14.0 

20.0 

18345-0641 

9.5 

1 

21.1 

20.9 

103 

40 

143 

1841 

24.0 

1.84 

15.0 

9.5 

12.0 

13.1 

18470-0044 

8.2 

1 

21.1 

20.6 

86 

16 

102 

1050 

11.0 

1.51 

14.8 

6.9 

7.1 

6.0 

18566-t0408 

6.7 

1 

20.9 

20.7 

25 

7 

32 

540 

9.1 

1.30 

7.5 

4.3 

7.2 

10.0 

19012-t0536 

8.6 

1 

21.3 

21.0 

105 

29 

134 

1339 

13.0 

1.25 

12.2 

11.0 

11.0 

8.9 

19012-t0536 

4.6 

1 

21.3 

21.0 

30 

8 

38 

383 

3.8 

0.67 

13.0 

2.9 

2.9 

2.4 

19035-t0641 

2.2 

1 

20.8 

20.3 

2 

1 

3 

28 

0.3 

0.32 

2.9 

0.9 

1.0 

0.9 

19217-tl651 

10.5 

1 

21.2 

21.1 

38 

70 

108 

1961 

36.0 

1.43 

7.5 

14.0 

26.0 

38.8 

19266-tl745 

10.0 

1 

20.9 

20.4 

17 

18 

35 

311 

2.8 

1.70 

18.4 

1.9 

1.7 

1.3 

19266-tl745 

0.3 

1 

20.9 

20.4 

0.02 

0.02 

0.04 

0.3 

0.003 

0.05 

1.0 

0.03 

0.03 

0.02 

19410-t2336 

6.4 

1 

21.2 

21.5 

92 

1331 

1423 

31896 

710 

1.33 

5.9 

240 

540 

982.9 

19410-t2336 

2.1 

1 

21.2 

21.5 

10 

143 

153 

3434 

77.0 

0.44 

3.8 

40.0 

90.0 

165.3 

19410-t2336 

6.4 

2 

21.2 

21.4 

96 

77 

173 

3766 

82.0 

1.40 

6.2 

28.0 

61.0 

108.3 

19410-t2336 

2.1 

2 

21.2 

21.4 

10 

8 

18 

405 

00 

bo 

0.46 

4.0 

4.6 

10.0 

18.2 

19411-t2306 

5.8 

1 

21.2 

21.1 

19 

28 

46 

649 

9.0 

0.98 

6.9 

6.7 

9.4 

10.8 

19411-t2306 

2.9 

1 

21.2 

21.1 

5 

7 

12 

162 

2.3 

0.49 

6.9 

1.7 

2.4 

2.7 

20216-t4107 

1.7 

1 

21.2 

21.1 

3 

3 

6 

43 

0.3 

0.26 

3.5 

1.7 

1.3 

0.8 

20293-t3952 

2.0 

1 

21.1 

21.5 

2 

7 

9 

270 

7.8 

0.31 

1.1 

8.6 

25.0 

59.2 

20293-t3952 

1.3 

1 

21.1 

21.5 

1 

3 

4 

114 

3.3 

0.20 

1.3 

3.0 

8.7 

20.6 

20293-t3952 

2.0 

2 

20.9 

20.2 

1 

0.1 

1.1 

21 

0.6 

0.15 

0.5 

1.6 

4.2 

9.1 

20293-t3952 

1.3 

2 

20.9 

20.2 

0.3 

0.1 

0.4 

8 

0.2 

0.09 

0.6 

0.6 

1.5 

3.2 

20343-t4129 

1.4 

1 

21.1 

21.2 

1 

1 

2 

9 

0.04 

0.16 

3.8 

0.6 

0.2 

0.1 

20343-t4129 

1.4 

2 

20.7 

20.8 

0.2 

0.2 

0.4 

1 

0.007 

0.17 

4.2 

0.1 

0.04 

0.01 

22134-t5834 

2.6 

1 

20.7 

21.2 

2 

15 

17 

242 

3.4 

0.47 

3.5 

4.9 

6.9 

7.9 

22570-t5912 

5.1 

1 

21.0 

20.9 

30 

43 

73 

911 

11.0 

1.48 

12.1 

6.1 

7.5 

7.8 

23033-t5951 

3.5 

1 

21.0 

21.2 

9 

23 

32 

566 

10.0 

0.85 

5.0 

6.4 

11.0 

16.9 

23139-t5939 

4.8 

1 

21.3 

21.1 

41 

16 

57 

662 

7.7 

0.67 

5.4 

10.0 

12.0 

11.7 

23151-t5912 

5.7 

1 

21.4 

20.8 

13 

8 

21 

597 

18.0 

0.55 

2.0 

10.0 

29.0 

72.7 


used for the low-mass flows discuss ed by [Cabrit &: Bertout 
(1992|)]and [Bontemps et al. (19961) is very simila r to ours 
and is based on results of Cabrit &: Bertout (199C| ). Hence, 
we can usefully compare our data with those results. 

Fig. 1^ shows the relations between the mechanical force 
Fqo, the core mass M^ore and the bolometric luminosity 
L. Fco is derived from Fm (Table ||) by applying an incli¬ 
nation correction. Because we do not know the inclination 
for each source separately we apply a mean correction fac¬ 
tor of 3, which corresponds to a mean inclination angle of 
57.3° ( Bontemps et ah, 1996 ). The main result of this 
comparison is that well established energetic correlations 
for low-mass outflows show a continuity up to the high- 
ma ss regime. A similar correlation was already presented 
by Lada (Ifl'S^ . 


The low-mass data can be divided into class 0 and class 
I sources, and Fig. ^(b) suggests tentatively that the ener¬ 
getic parameters in the high-mass regime correlate better 
with class I sources. This might be interpreted as an indi¬ 
cation that maybe no proper class 0 stage exists for mas¬ 


sive stars. But clearly this proposition has to be checked 
for larger samples to draw Arm conclusions. 

Nevertheless, the continuity between the low-mass and 
the high-mass regime is strongly suggestive of an outflow 
formation mechanism which is present both in the low- 
mass as well as in the high-mass case. 


4.2. Mass entrainment and accretion rates 


Fig. I presents an update of the mass entrainment rate 
versus bolometric luminosity p lot presented initially by 
Shepherd fc Churchwell (1996b ). In addition to our data 
we also include the massive outflow database compiled by 
[Churchwell ( 2000 ^ ) . Although the spread in the literature 
data is larger than for our more homogenous sample the 
mean of both is similar. This updated plot shows that the 
relation fitted by Shepherd &: Churchwell (1996b| ) to the 
data of Cabrit fc Bertout (1992| ) is just an upper envelope 
in the high-mass regime. The data indicate that the mass 
entrainment rate for sources > 1000 L© does not depend 
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Figure 4. Outflow parameter correlations: (a) mechanical force from CO versus core mass from 1.2 mm dust emission, 
(b) mechanical force versus bolometric luminosity. The dots represent our sample without distance ambiguity, the 
open triangles sources at near kinematic distance and the filled triangles sources at the far kinematic distance. The 
pentagons show data from Bontemps et al. (1996) and the squares data from [Cabrit fc Bertout (1992 ). Open squares 
and pentagons represent class I sources and fllled symbols class 0 sources. 


strongly on the luminosity, and has an average Mout = 
dMout/dt of 7.5 X 10“^ M 0 yr“^ with a spread between 
10 “'^ M 0 yr“^ and 10 “^ M 0 yr“^ (these numerical values 
are calculated just for our sample). While part of the non¬ 
correlation between M^ut and L at high luminosities is 
certainly due to a variety of observational uncertainties, 
we believe that much of the dispersion is real. One must 
remember that while outflows in low-mass star formation 
regions can usually be identifled with a given proto-star, 
in the regions of high-mass star formation studied here, we 
presumably deal with a cluster of embedded objects. It is 
possible that in some of our regions, for example, the most 
luminous object is on the ZAMS and has started nuclear 
burning whereas the most massive outflow emanates from 
a less evolved object. Proving this however will require 
higher angular resolution and comparison with VLA data. 

More interesting than the mass entrainment rates 
are the accretion rates Mace of such massive objects. 
Regarding the flows as momentum driven the momenta 
of the observed outflow and the internal jet entraining the 
outflow should be conserved, if there is efficient mixing at 
the jet/molecular ga s interface assumin g no loss of mo¬ 
mentum to the ISM ( Richer et ah, 2000 ): 


2 


Pflow — Mout'Cout — A/jet'Cjet^ — Pjet- 


■w 


L [logio(L0)] 

Figures. Mass entrainment rates versus luminosities 
(from [Sridharan et al. 2002| ). The dots (no distance am¬ 
biguity), open and filled triangles (near and far kine¬ 
matic distance, respectively) show the data presented in 
this paper, the squ ares represent the data compile d from 
IChurchwell ( 2000 b|) and the asterisks the data from Cabrit 


The mean value nflow of the observed nmax (Table ^ is ~ 
15 km/s. Correcting this by a mean inclination angle of 57° 
(multiply by 1/cos(57°)) results in a molecular outflow 
velocity of approximately 28 km/s. This is still a lower 
limit because we are noise-limited, and material at higher 
velocities is likely. 

.let velocities from high-mass proto-stars have only 


fc Bertout (1992 ). The dotted line shows the second or der 
polynomial fit to the data of Cabrit fc Bertout (1992| ) as 


proposed by Shepherd &: Churchwell (1996b|) . 


been determined in a few cases. Marti et al. (1998) ob¬ 


serve from proper motions jet velo cities ujet of a young 
massive star around 500 km/s, and Eisloffel et al. (2000) 
report that values between 500 km/s and 1000 km/s are 
typical. Based on this we assume in the following a mean 
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ratio between jet velocity and molecular outflow veloc¬ 
ity around 20 , which results in jet masses Mjet and mass 
loss rate of the jets Mjet about an order of magnitude 
below the rates of the entrained gas. Assuming further 
a ratio between mass loss rate of the jet and accretion 


rate of approximately 0.3 ( Tomisaka, 1998 ; Shu et al. 
199S(^e get a mean accretion rate of ~ 10 M 0 yr ^ for 
sources in the 10 ^ L© regime ; the distribution is shown in 
Fi gure |^. For the most luminous sources (~ 10^ L©) listed 
by Churchwell (2000b) the accretion rates get as high as 
10“^ MQyr“^. Such values are extremely high compared 
with typical accretion rates established for low-mass star 
formation (10“® — 10“® MQyr“^, Shu 1977). High mass 


accretion rates moreover ease the problem of overcoming 
radiation pressure and are requir ed in order to form mas - 
sive stars within a free-fall time ( Jijina fc Adams, 1996|) . 
It will be useful to examine the birth-li ne for such high ac¬ 
cretion rates in more det ail (see, e.g.. Stabler et al. 200C ; 
Norberg fc Maeder 2000 ). 



approximate accretion rates [Mq yr 


Figure 6. Histogram of the estimated accretion rates for 
our sample. 


Other interpretations are still possible, and we cannot 
rule out that the chaotic sources without clear outflow 
structure could be generated by interacting proto-stars. 
Collisions of proto-s tars are expected to be very energetic 
( Stabler et ah, 2000 ) causing either real explosions during 
stellar merger, or in less dramatic scenarios they could 
signiflcantly affect the rotation axis of the proto-stellar 
disks and thus cause strong precession of the outflows. 
But bipolar outflows are very difficult to explain in mas¬ 
sive star for mation scenarios based on the coalescenc e of 
proto-stars ( Bonnell et ah, 1998 ; Stabler et ah, 200^ ). In 
contrast, high accretion rates as found in this study give 
strong support to theories, which explain massive star for¬ 
mation mainly by accretion analogous to low-mass sce¬ 
narios albeit with accretion rates which are signiflcantly 
higher. The fact that highly collimated and massive flows 


are found (this work and studies with higher resolution 
by Beuther et al. (2002b ) and Wyrowski et al. in prep.) 
supports the proposition that the physical processes in the 
low- as well as the high-mass regime are similar. 


4.3. Outflow mass versus core mass 

We And a good correlation between the outflow mass Mout 
and the core mass Mcore derived from the 1.2 mm dust 


emission ( [Beuther et ah, 2002a|) , which is presented in Fig. 
0 (left). Our best fit over more than 3 orders of magnitude 
in core mass is Mout ~ 0.1 • M^ofe- The ratio Mout/Mcore 
is rather constant with an average of 0.04 for the sources 
(Figure right) and a spread of less than one order of 
magnitude. Thus, approximately 4% of the core gas is en¬ 
trained in the molecular outflow. There is a trend in our 
data that the ratio Mout/Mcore decreases with rising core 
mass Mcore, but the statistics at the upper mass end are 
too low for a definite conclusion on this point. 

From the correlation between outflow and core mass, 
we e stimate the accretion efficiency /acc- As discussed in 
§ 3.1 the derived parameters include contributions from 
several flows, but they are most likely do minated by one 
massive flow (see also the discussion in Tan & McKee 
2002 ), and we neglect effects due to multiple sources in 
the following. We define /acc 


/acc — 


Ma, 


M^core/^ff 


( 11 ) 


with the accretion rate Macc onto the star powering the 
massive flow observed by us, the free-fall time tg and the 
core mass Mcore- Furthermore, the ratio between the mass 
loss rate of the jet and the accretion rate is 


/r = 


Mi 


jet 


Ma. 


( 12 ) 


Multiplying eqs. ^ and using t rv ts (Fig. ||) we get: 


/r/acc = 


Mi 


jet 


Mjet _ Mjet Mout 


Mcore/^fF 


(13) 


The further assumption of momentum conservation (see 
eq. 1 ^ results in: 

'Cout M^out 


/r/acc = 


Ujet Me, 


(14) 


As outlined in §4^ we assume Uout/ujet to approximately 
1/20. Using additionally the average ratio between Mout 
and Mcore of 0.04 Eq. 14 reads: 


/r/acc = ^ X 0.04 ~ 2 X 10 


-3 


(15) 


The value 2 x 10 ^ in Eq. E has to be taken with cau¬ 
tion because the errors in core and outflow mass are both 
roughly a factor 2 , and it is not clear if the assumption 
of momentum conservation is correct or if some momen¬ 
tum is lost to the ISM exterior to the core observed at 
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Figure 7. The left plot shows the molecular outflow mass versus core mass ([Beuther et ah, 2Q02a| ) , whereas the right 
plot presents the ratio of both quantities versus the core mass. The circles show sources with no distance ambiguities, 
the filled and open triangles represent sources with far and near distances, respectively. 


1.2 mm. But it is interesting that the product /r/acc does 
not change significantly over many orders of magnitude 
in core and outflow mass. This implies that the ratio be¬ 
tween the ejected jet mass Mjet and the core mass Mcore 
is roughly constant for all outflows and cores (see eq. [T^ . 

4.3.1. Empirical estimates of /r and /acc 


Kicher et al. (2U(JU) describe an observational approach to 


estima te fr. They consider two options: 


(a) the luminosity is accretion dominated 

j. Fco ^'kep 
Jr — ^ Ukep 

L Vjet 


Hence, the derived parameters should be taken very cau¬ 
tiously, but nevertheless, Eqs. |^ & |l^ give approximate 
values for fr directly from the observations. Additionally, 
we can estimate from Eq. |T^ the accretion efficiency /acc- 
Table |3| gives the median values for fr and /acc taken 
via approach (a) and (b). Both scenarios - the accretion 
dominated and the ZAMS dominated - yield fr of order 
0 .1, which is smaller though not greatly (give n the un¬ 
certainties) than theoretical estimates of 0.3 by 


(1998|) or Shu et al. (1999 ). Richer et al. (2000 ) And mean 


Tomisaks 


(16) 


values fr around 0.3 in the accretion dominated case and 
values about an order of magnitude lower in the ZAMS 
case, which is consistent with our results. 


(b) the luminosity is mainly due to ZAMS stars 


fr — 


PCO 


IVL 


‘^kep 


(17) 


with Vkep the keplerian speed at the proto-stellar surface 
(equal to y/GM^/R^), which we assume to be the same 
as the jet velocity; Fco and pco are the mechanical lu¬ 
minosity and momentum from Table ^ multiplied by the 
above discussed statistical inclination correction factor 3, 
and M* and i?* are the stellar masses and radii for the 


ZAMS ( Lang, 1992 ). As luminosities, we use the bolomet- 
ric luminosities presented in Table |^. 


Th(i errorbars in both estimates are large (see also 


Richer et al. 2000 ). Masses and radii of the massive proto¬ 
stars are poorly constrained, and the momenta and me¬ 
chanical luminosities are uncertain to about an order of 


magnitude as well (Cabrit & Bertout, 1992). Thus, the es¬ 
timated errors in Eq. |16| fc |17| are larger than one order of 
magnitude. Additionally, the spread in the derived values 
of fr is about an order of magnitude for both approaches. 


Table 3. Median fr and face- (a) luminosity is accretion 
dominated; (b) luminosity is ZAMS dominated 

_ fr _/acc 

(a) 0.17 0.01 

(b) 0.07 0.03 

Rephrasing Eq. results in the accretion rate being 
a rather linear function of the core mass with the accre¬ 
tion efliciency and the free-fall timescale (^ 10^ yr) being 
approximately constant over the whole mass range. In the 
accretion do minated case (/acc ~ 0.01, see also ^ridharan 
et al. (200^ )) we And: 

Mace = 4^Mcore = 1 X 10“^ X Mcore [Meyr-^]. (18) 

Using additionally M* = Macc x t for the mass of the 
forming star, we get: 

= /accMcore 


(19) 
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As already stressed, the possible numerical uncertainties 
are large (about an order of magnitude), but equations 
111 & [^result, e.g., for a 5000 M© core, in an accretion 
rate of ~ 5 x 10“'* MQ/yr and a stellar mass of 50 M© 
for the most massive object of the evolving cluster. These 
are plausible values, and the accretion rate corresponds 
well to the range of accretion rates derived in §4.2. Such 
accretion rates are sufficiently high to overcome the radi¬ 
ation pressure and form massive stars (see, e.g., Jijina & 
Adams 1996). Additionally, relation agrees reasonably 
well with estimates of the most massive star of a cluster as¬ 
suming a star formati on efficiency of 30% (pLada, 1993 ) and 
a Miller-Scalo IMF (Miller & Scalo, 1979). We conclude 
that one can explain the formation of stars of all masses 
by similar accretion based scenarios with accretion rates 
roug hly proportional to final masses (Norberg fc Maeder, 
200C; i ionnell et ah, 2001 ; Tan &: McKee, 2002 ). 


5. Sunmary 


CO mapping of 26 sources showing line wing emission 
in previous single pointing observations reveal a high de¬ 
gree of bipolar outflow morphology with 21 sources being 
resolved into bipolar structures. This supports the idea 
that bipolar outflows are not only associated with low- 
mass star formation, but are ubiquitous phenomena for all 
masses. The size of the flows is on the pc-scale, and the 
data reveal collimation factors on the average > 2 , which 
is similar to results for low-mass flows but differs from pre¬ 
vious high-mass estimates in the literature (Richer et al. 


200C). We conclude that “Orion-type” flows of low col¬ 


limation factor are either rare in general or confined to 
very massive proto-stars (see however the discussion by 
Reipurth & Bally 2001). 

The dynamical timescales of the outflows correspond 
well to the free-fall timescales of the associated cores. This 
suggests that for high-mass flows, the CO timescale may 
often be a good estimate of the flow age. 

Derived flow parameters show that we are really deal¬ 
ing with massive and energetic sources, and compar¬ 
ing those quantities with equivalent studies in the low- 
mass regime shows that previously derived correlations 
for low-mass proto-stars have counterparts in the high- 
mass regime. We interpret that as support for star for¬ 
mation scenarios, which predict similar physical accretion 
processes for all masses with significantly increasing ac¬ 
cretion rates Macc in high-mass cores. 

We show that mass entrainment rates increase with 
core luminosity, and most of the massive sources show 
mass entrainment rates between ~ 10 “* M 0 /yr and ~ 
10“^ M 0 /yr with a few cases of higher values. In mo¬ 
mentum driven flows, these mass entrainment rates should 
correspond to mean accretion rates of order 10 “'* M 0 yr“* 
for sources of bolometric luminosity ^ 10^ L©. This is in 
good agreement with recent massive star formation sce¬ 
narios, which predict higher accretion rates for the most 
massive stars (Norberg & Maeder, 2000; Bonnell et al. 
2^; |Tan fc McKee, 200^ 


We find a tight correlation between the outflow mass 
and the core mass, which holds over more than three or¬ 
ders of magnitude in core mass. The ratio of both quan¬ 
tities is rather constant around 0.04. This correlation in¬ 
dicates that the product of the accretion efficiency and 
the ratio between the mass entrainment rate and the ac¬ 
cretion rate, which equals the ratio between jet and core 
mass /acc/r = Afjet/Afcore, is roughly constant during star 
formation of all masses. This makes the accretion rate (and 
thus the mass of the most massive star of the cluster) to 
first order a linear function of the core mass. 

Estimates of /r and /acc are around 0.2 and 0.01, re¬ 
spectively. In spite of large uncertainties, the results are 
consist ent with current jet entrainment flow formation sce¬ 
narios ( Richer et ah, 200C ). 

Those results support star formation theories in the 
high-mass regime which are based on similar principles as 
those for low-mass star formation (e.g., Norberg & Maedei 
_200C I lijina &: Adams 1996^ Wolfire & Cassinelli 1987; Tan 


fc McKee 20021). The intrinsic parameters are more ener¬ 


getic, but the physical processes are likely to be similar. 
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